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ABSTRACT 

We compare Ha, ultraviolet (UV) and infrared (IR) indicators of star formation rate (SFR) for 
a well-defined sample of z — 0.84 emission-line galaxies from the High-z Emission Line Sur- 
vey (HiZELS). Using emission-line, optical, IR, radio and X-ray diagnostics, we estimate that 
5-11 per cent of Ha emitters at this redshift are active galactic nuclei. We detect 35 per cent 
of the Ha emitters individually at 24 jim, and stack the locations of star-forming emitters on 
deep 24-jj.m Spitzer Space Telescope images in order to calculate the typical SFRs of our Ha- 
emitting galaxies. These are compared to the observed Ha line fluxes in order to estimate the 
extinction at z — 0.84, and we find a significant increase in dust extinction for galaxies with 
higher SFRs. We demonstrate that the relationship between SFR and extinction found in the 
local Universe is also suitable for our high-redshift galaxies, and attribute the overall increase 
in the typical dust extinction for z = 0.84 galaxies to an increase in the average SFR, rather 
than to a change in dust properties at higher redshift. We calculate the UV extinction, and find 
a similar dependence on SFR to the Ha results, but no evidence for a 2175 A UV bump in 
the dust attenuation law for high-redshift star-forming galaxies. By comparing Ha and UV 
indicators, we calculate the conversion between the dust attenuation of nebular and stellar 
radiation, 7, and show that 7 = 0.50 ±0.14. The extinction / SFR relationship is shown to be 
applicable to galaxies with a range of morphologies and bulge-to-disk ratios, to both merging 
and non-merging galaxies, and to galaxies within high- and low-density environments, imply- 
ing that it is a fundamental property of star-forming regions. In order to allow future studies 
to easily correct for a SFR-dependent amount of dust extinction, we present an equation to 
predict the extinction of a galaxy, based solely upon its observed Ha luminosity, and use this 
to recalculate the Ha luminosity function and star formation rate density at z = 0.84. 

Key words: infrared: galaxies, galaxies: evolution, galaxies: high-redshift, galaxies: ISM 



1 INTRODUCTION 

There are many methods of estimating the global star formation 
rate (SFR) of a galaxy, from measurements of the ultra-violet 
(UV) emission coming directly from visible young massive stars, 
to the bolometric infrared (IR) luminosity resulting from dust- 
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reprocessed stellar emission (see lKennicuttlll998l for a review of 
various star formation indicators). All SFR indicators suffer from 
different biases and uncertainties, and while SFRs estimated in 
different ways in the local Universe are now in general agree- 
ment, there is much great er scatter in the results at z ~ 1 (see, e.g. 
iHopkins & Beacomll2006h . 

Locally one of the best measures of the current rate of star for- 
mation within a galaxy comes from observations of the Ha transi- 
tion of atomic hydrogen. The luminosity of this recombination line 
scales directly with the rate of ionising flux from young, massive 
stars and gives an excellent indication of the presence of ongoing 
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star formation. However, its usefulness as a quantitative SFR in- 
dicator is reduced by the effect of dust obscuration; observed Ha 
line fluxes need to be corrected for extinction before the intrinsic 
Ha luminosity can be recovered, and the SFR calculated. 

The initial method used to estimate Ha extinction came from 
a comparison of thermal radio emission, which is unaffected by 
dust attenuatio n, and the observe d Ha + [N II] line emission within 
local galaxies lKenn icut3[l983n . Disentangling thermal and n on- 
thermal radio emission is non-trivial (see, e.g. ICondorjfl992l for 
more details), as thermal radio emission is only a small fraction of 
the total radio lumin osity of a galaxy (< 10 per cent at 1.4 GHz; 
ICondon & Yin|[l99oT) , and so this method is not often used in prac- 
tise. 

A much more common method to estimate extinction in in- 
dividual galaxies comes from the 'Balmer decrement', compar- 
ing the observed value of the Ha/H/3 line flux ratio to the in- 
trinsic value when no dust is present. This allows a calculation of 
the reddening between the wavelengths of the Ha a nd H/3 lines 
and, with a choice of dust attenuation model (e.g. Calze tti et al.l 
2000), permit s an estimate of the attenuation at a g iven wave- 
length (see e.g. lMoustakas. Kennicutt & Tremontil 20061) . However, 
this technique requires sensitive spectroscopy to obtain the line 
flux ratio with a good signal-to-noise ratio (SNR), which is time- 
consuming for large samples of galaxies at moderate redshifts. It 
is therefore more usual for studies with limited available photom- 
etry to simply assume that a single extinction can be applied to 
all galaxies. The canonical amount of Ha extinction for a typi- 
cal g a laxy is usually taken to be Ah« ~ 1 mag ( e .g. iKennicutt 
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of extinction in the local Universe. 

At higher redshift, the appropriate Ha extinction for a typ- 
ical galaxy is not well determined. Significant numbers of ex- 
tremely dusty galaxie s have been discovered at high redshift (e.g. 
Chap man etail 2005). and the large dust content implies a greater 
Ha attenuati on. The global rate of st ar formation increases out 
to z ~ 1 (e.g. iHopkins & Beacomll2006h . and galaxies with higher 
SFR are thought to undergo greater ext inction iHopkins et alj200ll : 
ISullivan etld1l200ll ; iBerta et al.ll2003l) . also implying that the typ- 
ical H q attenuation will be greater at z ~ 1 than in the local Uni- 
verse. IVillar et al]fe008h have estimated that the typical Ha extinc- 
tion of galaxies within their z = 0.84 sample is ~ 1.5 mag. 

The total IR luminosity of a galaxy is a bolometric measure- 
ment of the thermal emission from dust grains, which have been 
heated by t he UV and optic al radiation produced by young massive 
stars (e.g. Kennicutt 1998). For star-forming galaxies with a rea- 
sonable dust content, this emission gives a direct measurement of 
the total stellar power output. However, since longer wavelength IR 
emission will more closely trace the cool dust component within a 
galaxy, which predominantly results from heating by an old stellar 
component rather than from recent star formation, it could be ar- 
gued that mid-IR emission mi ght be a better SFR indicator than the 
bolometric IR luminosity (e.g. lDopita et al.l2002h . 

iRieke et al] d2009h have demonstrated that the flux density 
measured in the Spitzer Space Telescope 24- flm band can be used 
as a good proxy for SFR for galaxies at z = 0, and, after applying 
suitable corrections, has a comparable accuracy to estimating SFRs 
from the extinction-corrected Paa luminosity or total IR luminos- 
ity for normal galaxies in the redshift range < z < 2. Observa- 
tions at 24 flm have the additional advantage over longer wave- 
length IR observations in that they are more sensitive (thus allow- 



ing fainter, more distant galaxies to be studied), and have a higher 
resolution (~ 6 arcsec, making it relatively straightforward to iden- 
tify counterparts at other wavelengths; the 70- and 160-flm Spitzer 
bands only have resolutions of ~ 20 and ~ 40 arcsec respectively). 
These disadvantages will be overcome to some extent with Her- 
schel, which should be able to make sensitive ~ 6 arcsec resolution 
images of galaxies at 70 flm, close to the peak in luminosity for the 
spectral energy distribution (SED) of a typical galaxy at z ~ 0. 

In this work, we compare Ha, UV and IR SFR indicators in 
order to estimate the extinction for a well-defined sample of galax- 
ies at z = 0.8 4, using data from the High-z Emissio n Line Sur- 
vey (HiZELS: iGeach et alj|2008l ; ISobral et alj|2009allbl) . HiZELS 
is a panoramic extragalactic survey making use of the wide area 
coverage of the Wide Field Camera on the 3.8-m UK Infrared 
Telescope (UKIRT). Narrow-band observ ations of the Cosmo - 
logical Evolution Survey field (COSMOS: [Scoville et al]|2007al) . 
and the Subaru-XMM -UKIDSS Ultra Deep Survey field (UDS; 
lLawrence et alj 12007) have been taken in the J, H and K bands 
in order to target Ha emitters at z = 0.84, z = 1 .47 and z = 2.23 re- 
spectively; this paper concentrates on the properties of the z = 0.84 
Hq sam ple, the creation of which is described in full in Sobr al et"aD 



In Section [2] we present details of the HiZELS sample selec- 
tion, and the available multi-wavelength data in the COSMOS and 
UDS fields. In Section[5] we describe various methods which we use 
to identify active galactic nuclei (AGN) contaminants, and quantify 
the fraction of Ha-emitting sources which may be AGN at z = 0. 84. 
We describe our method for stacking galaxies and measuring typ- 
ical 24- flm flux densities in Section [4] and convert our stacking 
results into SFRs in order to estimate the Ha extinction that is 
required to make this indicator agree with the IR-based SFRs. In 
Section[5]we discuss our results, investigate the effects that a vari- 
able extinction has on the Ha luminosity function at z = 0.84, and 
present a method for correcting observed Ha luminosities for a 
SFR-dependent amount of extinction. Section[6]reports our conclu- 
sions. 

Throughout this work we assume a concordance cosmology of 
SlM = 0-3, Q. A = 0.7 and H p = 70 km s~' Mpc ~'. All magnitudes 
are given in the AB system dOke & Gunnlll983h . 



2 HIZELS SAMPLE SELECTION 

We give a brief overview of the HiZELS sample selec t ion pro - 
cedure below; full details can be found in ISobral et all d2009al) . 
Emission-line candidates were selected on the basis of a flux excess 
in a narrow /-band filter centred on A = 1.21 1 flm, compared with 
deep ./-band continuum observations. Spurious sources (mainly ar- 
tifacts due to bright stars) were cleaned from the sample, and spec- 
troscopic and photometric redshifts were used to identify and re- 
move sources found to be a different line emitter at another redshift 
(such as Hj8 or [Olll] at z ~ 1.5). 

The Ha line flux for each source was calculated from the 
measured flu x densities i n the na rrow and broad-band filters, as 
described in ISobral et alj {2009a). No aperture corrections were 
applied to the Ha flux, since the correction factor is expected to 
be minimal. Contamination from the nearby [Nil] lines at 6548 
and 6583 A was corrected for, using the relationship between 
the flux ratio Ft^m/F^ and the total measured equivalent width 
EW(Ha+[Nll]) from lVillar etaD d2008h . with the median correc- 
tion factor being 25 per cent. 

The final HiZELS Ha sample consisted of 477 Ha-emitting 
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sources in COSMOS and 266 in UDS, over a total area of 
1.30 deg 2 , at a redshift of z = 0.845 ± 0.021, down to an average 
observed 3(7 Ha line flux limit of 8 x 10~ 17 erg s -1 cm~ 2 . 



2. 1 The HiZELS COSMOS sample 

Ground-based observations of the COSMOS field have been 
made at a variety of UV, optical and near-IR wavelengths, and 
ICapak et all 420071) have produced a point-spread function (PSF)- 
matched catalogue for all sources with /-band magnitude < 25. 
This catalogue contains data from observations with the Canada 
France Hawaii Telescope (CFHT; u* and bandfl), the Sl oan Dig- 
ital Sky Survey (SDSS ugriz bands; lAbazaiian et alj|2004l) and the 
Subaru telescope (Bj, Vj, g + , r + , i + , z + and NB816 bands), along 
with A" s -band data taken from a combination of observations from 
the Cerro Tololo International Observatory and the Kitt Peak Na- 
tional Observa tory. Further detail s on the catalogue construction 
can be found in lCapak et ail d2007T) . All of the COSMOS Ha emit- 
ters have counterparts in this PSF-matched catalogue. 

Mid-IR and far-IR observations of COSMOS hav e been taken 
as par t of the Spfeer-COSM OS survey (S-COSMOS; ISanders et ail 
l2007l ; lLe Floc'h et al.ll2009T) . 464 of the Ha emitters (97 per cent) 
have a counterpart in Infrared Array Camera (IRAC) observations 
of the field, in at least one of the 3.6-, 4.5-, 5.8- and 8-/lm bands. 
Several of our AGN classification diagnostics rely upon the avail- 
ability of IRAC information (see Section [3}; for this reason, we 
reject the 13 sources without IRAC detections from our study. To 
determine whether there is any potential bias from this exclusion, 
we repeated all tests in Section|4]while retaining these sources - no 
significant differences were seen in our results. 

There have been several Spitzer 24-/xm surveys of the COS- 
MOS field, and we have used the most recent of these, the Multi- 
band Imaging Photometer for Spitzer (MIPS) Cycle 3 observa- 
tionfl in this work. The public source catalogue is cut at a flux 
density limit of 0. 15 mjy (~ 7(7), and 158 Ha emitters had a coun- 
terpart in this catalogue, where we consider that an Ha emitter has 
a counterpart if the positions of the 24-flm and optical source cen- 
tres are within 3 arcsec. We estimate the percentage of Ha emitters 
which may have an incorrect 24-flm association by shifting their 
positions by +10 arcmin in declination and recalculating the num- 
ber of matches - nine spurious associations were made, which im- 
plies an overall percentage of Ha emitters with individual 24-/im 
counterparts of 34 ± 2 per cent. Throughout this work, all matching 
percentages have been corrected for false matches in this way. 

1 .4-GHz radio observations of the COSM OS field have been 
taken with the Very L arge Array (VLA; ISchinnerer et alj 1 20041 
l2007l : lBondietal.ll2008h . The Ha emitters are located within the 
VLA-COSMOS Deep Project area, which has a typical noise 
level of between 10 and 40 /ijy beam -1 , and 21 Ha emitters 
(4.5 ±0.9 per cent) have 1.4-GHz counterparts within 3 arcsec. 
The full HiZELS reg ion of the COSMOS fi eld has been observed 
with XMM-Newton dCappelluti et al.l [2009). and the majority of 
the region has been observed with Chandra jElvis et al.ll2009f) - 
7 (1.5 ± 0.4 per cent) and 8 (1.7 ± 0.2 per cent) Ha emitters have 
X-ray counterparts within 3 arcsec respectively, with 5 Ha emitters 
having counterparts from both surveys. 



1 Throughout this work we retain the notation of ICapak et alj b007t) to 
differentiate between observations in different filter systems. 

2 available via http://www.ifa.hawaii.edu/~ilbert/S-Cosmos/ 



2.2 The HiZELS UDS sample 

The UDS field is the locat ion of a very deep U KIRT near-IR survey 
in the J, H and K bands JWarren et ai1l2007l) . There is additional 
deep imaging data av ailable from Subaru (Bj, Vj, r + , v, z + bands; 
iFurusawa et al.|[2008h . and «*-band data from CFHT. 

A deep Spitzer legacy survey of the UDS field has been carried 
out with both IRAC and MIPS (SpUDS; P.I. J. Dunlop). Catalogue 
production is still ongoing, but mosaics of the SpUDS survey have 
been delivered to the communit^ We use a preliminary version 
of the IRAC catalogue (M. Cirasuolo et al., in preparation), which 
excludes regions of the UDS image near to bright stars, image ar- 
tifacts or the mosaic edges - 234 of the Ha emitters are in regions 
covered by this catalogue. We exclude a further four sources from 
our analysis because they are either outside or on the edge of the 
SpUDS 24-flm image, leaving 230 sources, all of which had optical 
/ IR counterparts. Photometry for all bands was measured within a 
1.8 arcsec diameter aperture, equivalent to roughly 2x the size of 
point spread function for the optical and near-IR images. 

We make use of the finely sampled 24-flm SpUDS image of 
the field (1.245 arcsec/pixel). This image has asteroid streaks visi- 
ble; we inspected the locations of the Ha emitters to confirm that 
these imaging artifacts would not affect our analysis. 84 of the UDS 
emitters (37 ± 3 per cent) had counterparts in the preliminary 24- 
flm source catalogue, which we have cut at 0. 1 mjy (~ 6(7), using 
the same matching criteria as for the COSMOS field. 

Radio observations of the UDS field have been taken at 
1.4 GHz l lSimpson et all f2006). and three of the Ha emitters 
(1.3 ±0.4 per cent) have 1.4-GHz counterparts above 100 ^iJy. 
XMM-Newton observations of the fiel d have been taken as part of 
the larger Subaru/XMM Deep Survey dUeda et al.l ! 2008). and four 
of the Ha emitters (1 .7 ± 0.4 per cent) have X-ray counterparts. 



3 IDENTIFICATION OF AGN CONTAMINANTS 

The HiZELS selection process identifies Ha emitters within the 
redshift range of z = 0.845 ±0.021. While most of these sources 
are likely to be powered by star formation, there will be some ob- 
jects where AGN activity dominates the energetics - SFRs esti- 
mated from these sources will not be accurate, and these contami- 
nants need removing from the sample. We adopt a number of meth- 
ods for classifying sources as 'possible' or 'likely' AGN, with the 
number of sources satisfying each of the criteria listed in Table Q] 

3.1 SED template fitting 

The matched optical and IR photometry available for each Ha 
emitter was described in Section [2] This photometry was blue- 
shifted to the source re st-frame and fitted to a series of tem plate 
SEDs (shown in Fig. [Tl bale et al.ll200ll ; IPolletta et alj|2007l) . us- 
ing a standard x 2 minimization procedure. Only the data with the 
highest SNR was used out of multiple observations at similar wave- 
lengths, such as the SDSS u and CFHT u* bands. The SEDs con- 
sisted of star-forming galaxies (spiral, elliptical and starburst tem- 
plates), Seyferts, and AGN templates (Type-1, Type-2, and two 
composite starburst/ AGN templates). Sources which were best fit- 
ted by an AGN SED were classified as 'possible' AGN candidates 
- we do not classify sources as 'likely' AGN through this automatic 
fitting procedure, as the range of SED templates is limited to a few 

3 available via http://ssc.spitzer.caltech.edu/legacy/spudshistory.html 
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Rest-frame wavelength (//m) for z = 0.84 

Figure 1. The te mpl ate SEDs used to fi t the Ha emitters, taken from 
iDale et alj <200ll) and lPolletta et alj 120071) and shown in the galaxy rest- 
frame, arbitrarily normalised to a flux density of 1 at 5500 A. The vertical 
dashed lines indicate the effective wavelengths of the available photometry 
for sources at z = 0.84, with the narrow/broad-band J selection shown by 
the thicker grey line. Note the difference in AGN and star-forming galaxy 
mid-IR SED shapes. 



per source type, and may not fully represent the HiZELS popula- 
tion. Fig.[T]shows that while there is a large amount of photometric 
information available for the Ha emitters, there is very little dif- 
ference between the colours of the Seyfert and Starburst templates 
where photometry is available - for this reason, sources that were 
best-fitted by Seyfert templates were excluded from this automatic 
classification. 

The SEDs of all sources were inspected to confirm the results 
of the automatic fitting, and to identify any further objects which 
appeared to be potential AGN, but had not been classified as such. 
Four further sources were manually classified as possible AGN; 
these typically had a similar shape to the AGN templates but had 
an offset in the IR to optical flux densities, and for this reason were 
not automatically classified as AGN. 

3.2 Colour-colour diagnostics 

We use a number of colour-colour diagnostics to identify further 
AGN contaminants. For each, we classified sources located within 
the selection region as 'likely' AGN, while those located within 
the region, but with errors on their flux densities sufficiently large 
enough to potentially place them outside the region, were classified 
as 'possible' AGN. 

We used the mid-IR selection criteria of Istern et alj {2005) to 
define a region of IRAC colour-colour space that principally con- 
tains AGN, which is shown on Fig.[2ja), along with the positions 
of the Ha emitters with detections in all four IRAC bands. The ob- 
jects which have been flagged as potential AGN contaminants lie 
well away from t he re mainder of the population. 

Stern et al. I fe005t) warn that their criteria may preferentially 
omit AGN at z — 0.8, where the difference in mid-IR colours be- 
tween star-forming galaxies and AGN becomes less distinct. Fig.Q] 
shows that the composite black-body spectrum of galaxies powered 
by star formation peaks at around 1.6 flm, and falls at longer wave- 
lengths. Warm dust located near to an AGN increases the luminos- 
ity of the mid-IR continuum, and Fig. Q] indicates that the AGN 
templates continue to rise with wavelength throughout the mid-IR. 
For this reason we add a further AG N classification for sources that 
are outside of the lStern et al.1 J20051) region, but have a rising SED 
between the 3.6- and 4.5- flm bands (which probe rest-frame 1.9- to 
2.4- flm dust emission), which is also shown on Fig.|2ja). 

It is difficult to conclusively discriminate between Seyfert and 



Table 1. The number of Ha emitters classified as 'possible' or 'likely' AGN 
under each of the criteria described in Section fj] Several sources fulfilled 
more than one criterion, so the total number of sources classified as AGN 
is less than the sum of the individual classifications. Sources with two or 
more 'possible' classifications were considered to be 'likely' AGN. Values 
in parentheses represent the number of sources which were only flagged as 
AGN under a single criterion. 



Method 


Type 


COSMOS 


UDS 


Total 


Automatic SED fit 


Possible 


17 (10) 


20 (10) 


37 (20) 


Manual SED fit 


Possible 


4(4) 


(0) 


4(4) 


Stern wedge 


Possible 


3(2) 


17(8) 


20 (10) 


Stem wedge 


Likely 


7(0) 


2(0) 


9(0) 


[3.6] - [4.5] colour 


Possible 


1(1) 


2(1) 


3(2) 


[3.6] - [4.5] colour 


Likely 


4(3) 


0(0) 


4(3) 


Mid- to far-IR slope 


Possible 


3(3) 


2(0) 


5(3) 


Mid- to far-IR slope 


Likely 


9(1) 


5(0) 


14(1) 


X-ray counterpart 


Likely 


10(3) 


4(0) 


14(3) 


Low value of 514 


Likely 


1 (0) 


3(3) 


4(3) 


Emission-line ratio 


Likely 


Ss 3 (2)" 




> 3 (2)" 


Ha emitters 




464 


230 


694 


Possible AGN 




20 


19 


39 


Likely AGN 




22 


15 


37 


AGN contamination 




5-9% 


7-15% 


5-11% 


Sources retained for analysis 


422 


196 


618 



" This value is a lower limit, as only 26 COSMOS sources (6 per cent) had 
spectroscopy with sufficient SNR for this method. 

starburst galaxies, given the available photometry - the optical SED 
shape of the two source types can be very similar, with the main 
difference being the increased amount of 24- flm emission for star- 
bursts. We use the shape of the mid-IR to far-IR SED to define a re- 
gion R(5.S, 8, 24) that contains Seyferts, from inspection of Fig.Q] 

#(5.8,8,24) = ([8] - [24] < 2 and [5.8] - [8] > -0.1). (1) 

The first criterion distinguishes between Seyferts and starbursts, 
and the second criterion prevents more quiescent source types such 
as spirals and elliptical galaxies from being flagged. Fig.[2lb) shows 
the selection region for these criteria, which also includes the other 
AGN SEDs shown on Fig.0 

3.3 X-ray detections 

The COSM OS field has bee n observed w ith both XMM-Newton and 
Chandra dUeda et alJkOQSt lElvis et al.ll2009h. and the UDS field 
has been observed with XMM-Newton (Cappellutiet al. 2009). The 
catalogue flux limits are equivalent to luminosity limits in the 2- 
10 keV energy band of 3 x 10 43 erg s _1 , 2.5 x 10 42 erg s _1 and 
1 x 10 43 erg s _1 respectively, for sources at z = 0.845. Purely star- 
forming galaxies can emit X-rays, but locally are not observed 
to have luminosi ties exceeding Lo 5_g j^-v ~ 3 x 10 42 erg s -1 
jBauer etal ■120041) - all of the Ha emitters with an X-ray detection 
have luminosities significantly above this value, and are classified 
as 'likely' AGN. 

3.4 Radio-loud sources 

We use the IR / radio correlation (e.g. lAppleton et al.1 l2004t) to 
distinguish between radio emission resulting from AGN and star- 
formation activity, for the Ha emitters with radio counterparts. 



Obscured star formation at z= 0.84 5 




Star-forming 
AGN 
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[5.8] - [8.0] (AB) [5.8] - [8.0] (AB) 

Figure 2. Infrared colour-colour plots showing the HiZELS Ha emitters (COSMOS galaxies with black dots, UDS galax ies with grey dots ; AGN candidates 
with black diagonal crosses and grey upright crosses respectively). AGN selection regions are indicated by the lines: (a) the Stern et al. (2005) criteria is shown 
with the solid line, and the additional 3.6- to 4.5-/im criterion is indicated by the dashed line; (b) the region R(5.8,8,24) is indicated by the dashed line. The 
colours of the SEDs shown on Fig.[T]are also plotted for reference, separated into star-forming galaxies (solid black squares), Seyferts (open triangles) and 
AGN (open circles). 



22/24 of these had catalogued 24-flm counterparts - visual inspec- 
tion of the MIPS images confirmed the lack of a counterpart for the 
remaining two sources, and we place a conservative upper limit of 
0.1 mjy on the 24-flm flux density of these sources. 

The logarithmic ratio of 24-flm and 610-MHz flux densities 
is 0.55 ±0.4 for star-forming g alaxies at z ~ .8 which are de- 
tected at both wavelengths (see iGarn et al ] 12001 for further de- 
tails); converting this to a 1.4-GHz value using equation (4) of 
IGarn & Alexander! d2009l) . and assuming a radio spectral index of 
0.8, we obtain a typical value of ^24 = log 10(^24 /Si .4) = 0.8±0.4 
for star-forming galaxies at z = 0.84. We find that indeed most of 
the HiZELS radio emitters do cluster around <5>24 ~ 0.7, but four 
lie well away from this, with ^24 < 0.1, and these are classified as 
'likely' AGN. 

3.5 Emission-line ratios 

Spectroscopic data is available for 20 per cent of the COS MOS 
Ha emitters from th e z-COSMOS survey dLillv et alj|2009h . and 
ISobral et al.1 fe009al) used the [O Il]3727/H/3 and [O III]5007/H/J 
emission-line ratios to separate AGN from star-forming galaxies. 
26 of the emitters had spectroscopy at sufficient SNR to carry out 
this co mparison; of t hese so urces, three (12 per cent) were classi- 
fied bv lSobraletal] {2009a) as likely AGN contaminants, and we 
retain that classification in this work. No spectroscopy is currently 
available for the UDS HiZELS sample, and as we do not have com- 
plete spectroscopy for the entire COSMOS sample, this value is a 
lower limit. 



as AGN through the IR colour-colour criteria; 10/464 (2. 1 per cent) 
of COSM OS sources and 19/2 30 (8.3 per cent) of UDS sources fell 
within the IStern et al.l (2005J) region. We have placed more con- 
servative errors on the UDS photometry as production of a band- 
merged catalogue is still ongoing - this leads to a higher fraction of 
UDS sources being 'possible' AGN than in COSMOS. 

We visually inspected the SEDs of these sources in order to 
investigate this discrepancy - the majority showed a clear power- 
law AGN shape, and were also flagged as AGN under at least one 
other criteria, making it very likely that the classifications are cor- 
rect. We believe that the bulk of the discrepancy between the two 
fields may be due to cosmic variance, although the slightly greater 
errors on the UDS photometry (which may incorrectly move a few 
sources into the AGN classification regions) may exaggerate this. 

By combining the samples, we obtain an overall AGN contam- 
ination rate of between 5 and 1 1 per cent. The majority of sources 
classified as 'possible' AGN contaminants were best-fitted by an 
AGN template, but did not satisfy any of the specific colour require- 
ments; for this reason, we believe that the true AGN contamination 
rate will be near to the upper end of our quoted range. The varia- 
tion in the number of sources detected by each technique suggests 
that (i) using a single criterion for starburst / AGN discrimination 
will underestimate the true AGN fraction, and (ii) a combination of 
colour-colour methods and full SED fitting should be used in order 
to maximise the success rate for AGN identification. 



3.6 The AGN contamination rate at z = 0.84 

Several of the Ha emitters were classified as AGN under more than 
one criteria, and sources with two or more 'possible' classifications 
were considered to be 'likely' AGN candidates. In practice the dis- 
tinction between 'possible' and 'likely' AGN candidates is only of 
interest in quantifying the uncertainty in the fraction of Ha emit- 
ters which may be AGN; we do not use either of these types of 
source in later analysis. The total 'possible' and 'likely' AGN in 
each field were 20 and 22 (COSMOS) and 19 and 15 (UDS). We 
find a significantly different AGN fraction between the two fields 
(5-9 per cent in COSMOS, and 7 - 15 per cent in UDS). The prin- 
cipal cause of this discrepancy is the number of sources classified 



The technique of source stacking (used in Section @ im- 
plicitly assumes that the sources being combined statistically are 
comparable to each other, but it is robust to the presence of a 
small number of outlier sources, such as AGN contaminants (see 
IGarn & Alexander 2009 ; , for further details). The results in the re- 
mainder of this work use only those sources that were not classi- 
fied as either likely or possible AGN, but we have repeated all tests 
while retaining the possible AGN in our sample, and found no sig- 
nificant difference in any of the results. The Ha and 24-flm flux 
density of the AGN contaminants are very weakly correlated (see 
Section l4~Tt . and we believe that it is unlikely that any remaining 
AGN contaminants should have a significant effect on the results 
presented in this work. 
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Figure 3. Variation in 24-ilm flux density with Ha flux, for the stacks of all Ha emitters which were not flagged as AGN. The solid diagonal line represents the 
best-fitting relationship between 24-ilm flux density and Ha flux, with a power-law slope of 0.83. Individual Ha emitters with 24-jim detections are indicated 
for comparison with the stacked data, separated by their AGN classification, and the COSMOS catalogue limit of 0.15 mJy (~ 7o") and UDS flux limit of 
0.1 mJy are indicated by the horizontal solid and dashed lines respectively. Note that the filter profile correction, described in Section l4~71 has been applied to 
the stacked Ha fluxes but not to the individual data po ints. The right-h and axis indicates the corresponding SFR for each 24-itm flux density, calculated using 
equation {5). The diagonal dotted lines represent the Kennicuul H998h relationship between Ha luminosity and SFR, modified by a factor of 0.66 in order to 
make the IMFs agree (see Section l4~2l . and assuming an extinction of 1, 2, 3 and 4 mag (bottom to top). 



4 ANALYSIS & RESULTS 

4.1 Calculation of stacked 24-fim flux densities 

It is necessary to know the locations of sources to within the size 
of a pixel (1.2 and 1.245 arcsec for the COSMOS and UDS 24- 
images respectively) in order to use source stacking. We have used 
source positions taken from the optical catalog ues, which are ex- 
pected to have an accuracy of about 0.2 arcsec dCapak et al]|2007l ; 
iFurusawa et al.|[2008t> . and we verified that there is no significant 
coordinate offset between the HiZELS data and their counterparts 
in the 24-flm images (median offsets in right ascension and decli- 
nation of -0.04 ±0.04 arcsec and 0.03 ±0.04 arcsec for COSMOS 
sources, and —0.05 ±0.08 arcsec and 0.30 ±0.09 arcsec for UDS 
sources). The small positional offsets were not corrected for in later 
analysis. 

The stacking technique used is very similar to that described 
in I Garn & A lexander (2009), and further details and validation of 
the procedure can be found in that work. Sources were binned by 
their observed Ha line flux, and a median Ha flux calculated for 
each bin. A typical 24- fim flux density was then calculated for each 
set of binned data in two different ways: 

(i) The 24-flm flux density of each Ha emitter was measured 
within a circular aperture of 3 arcsec radius, centred on the source 
location. While individual flux density measurements are affected 
by the local 24-flm noise level, the statistical properties of the 
distribution give a robust estimate of the typical flux density for 
sources within the bin. The median was used as the statistical esti- 
mator, to avoid being b iased by the combination of detected and un- 
detected sources (see lGarn & Ale xander 2009, for further details), 
as for many stacks, the median flux density will come from indi- 
vidual detections at 24 flm (see Fig. [3}- 



(ii) A cut-out 24-flm image was created for each source in a 
bin, centred on the known HiZELS locations. From these cutouts, 
a stacked image was created, with each pixel calculated from the 
median of the values of the corresponding pixels of the cutout im- 
ages. The flux density of the stacked image was then measured in 
the sa me way as for method (i). We experimented with the sugges- 
tion of lLonsdale et alj (2009) of placing images into a stack twice, 
once at 90° to the original orientation, in order to remove any sys- 
tematic instrumental effects, but found that this made a negligible 
difference to the stacking results. 

Measurements of the flux density obtained from random posi- 
tions indicated the presence of a slightly negative background level 
in both images (approximately —0.01 MJy sr~', equating to a flux 
within a 3 arcsec radius of —8 flly), which was subtracted from the 
flux density measurements to avoid biasing the stacking results. 

The measured 24-flm flux density was multiplied by an 
aperture correction factor, obtained through comparison of the 
measured and catalogued flux densities of 24-flm sources 
within the COSMOS fie ld, using a similar method to that in 
iGarn & AlexandeJ {2009). The calculated correction factor was 
3.41 ±0.12, and should also be appropriate for UDS - the agree- 
ment between the stacked results from the two fields (see Fig. [3} 
implies that we are not introducing a systematic effect by using the 
COSMOS-derived aperture correction for the UDS stacking. 

The narrow 7-band filter used to identify Ha emitters is not 
a perfect top-hat function, and this can have a measurable effect 
for sources with Ha emission at wavele ngths near the edge of the 
filter (as discussed in Sobr al et al . 2009a). We calculate the statisti- 
cal effec t of this by simulati ng a distribution of 10 5 sources based 
upon the lSobral et al.1 d2009al) Ha luminosity function. We spread 
them evenly in redshift over the range z = 0.81 — 0.87, and use the 



Obscured star formation at z = 0.84 7 



actual filter profile to recover the luminosity distribution. Through 
comparison of the median input and median recovered luminosity, 
a statistical correction to the observed luminosity of each source 
was obtained - these corrections were applied to the Ha data be- 
fore stacking took place. The filter profile correction has a moderate 
effect on the faintest emitters (up to a ~ 20 per cent increase in lu- 
minosity), and little effect for the brightest emitters, as these can 
typically only be observed to be bright if they are not detected in 
low transmission regions of the filter. Nevertheless, all of the trends 
discussed in this work remain present if no correction is applied. 

Error estimates for the Ha fluxes and 24-flm stacked flux 
densities indicate the 68 percent confi dence limits (+lo ~) obtained 
from the measured flux distribution fc}ottetal.ll200lh . The flux 
densities obtained using both stacking methods were consistent to 
within 3 per cent, which we consider to be an additional uncertainty 
introduced through the stacking process, and add in quadrature to 
our error estimate. Throughout this work, all stacking results have 
been taken using method (i). The dominant errors on our stacked 
results come from the uncertainty in estimating stacked flux den- 
sities, due to the small number of galaxies in each bin (typically 
between ~ 20 and ~ 160). 



At z = 0.845, the Rie ke et alj d2009h relationship is 



4.2 Variation in 24-jUm flux density or SFR with observed 
Ha flux 

Fig.|5]shows the stacked flux density measurements using the com- 
bined data from the COSMOS and UDS fields, with the results 
from each individual field also indicated for reference. While there 
is a strong correlation between the two stacked flux measurements, 
the relationship is significantly non-linear; the best-fitting power 
law has a slope of 0.83 ±0.05, and is indicated on Fig. [3] We 
have plotted individual Ha emitters with 24-flm detections on 
Fig. [3] separated by their AGN classification, in order to investi- 
gate whether there is any strong dependence of 24-flm flux density 
on Ha flux for the AGN sources with individual 24-flm detections 
- no such dependence is seen, with a Pearson product-moment cor- 
relation coefficient of p = 0.29, compared with p = 0.55 for the 
sources that are not classified as AGN. 

Riek e et al .(2009) present a method for determining SFRs di- 
rectly from Spitzer 24-flm flux density measurements. They assem- 
ble a series of template SEDs for galaxies with different luminosi- 
ties, and create average SED templates which span the total IR lu- 
minosity range of 5 x 10 9 — 10 13 Lq. From these bolometric lumi- 
nosities^ SFR i s calculated using a similar relationship to that in 
iKennicutl dl998l) . Each template is redshifted in steps of Az = 0.2, 
and convolved with the MIPS 24-flm instrumental response func- 
tion, in order to calculate the flux density which would be observed 
for the galaxy at a given redshift. The relationship between ob- 
served 24-flm flux density and SFR is parameterised at each dis- 
crete redshift, with the relationship at intermediate red shifts avail- 
able th rough interpolation of the fit coefficients. The Rie ke et ail 
d2009h parameterisation takes into account the variation in spectral 
shape (i.e. the K-correction) which is needed when converting from 
an observation-frame 24-flm flux density into a rest-frame 24-flm 
luminosity. Accurate knowledge of the K-correction is very impor- 
tant for 24-flm observations, which do not just probe the hot dust 
emission from the galaxy, but can also detect different pol ycyclic 
arom atic hydrocarbon features (PAH molecules; see e.g. iDrainel 
120031) at varying redshifts. 
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where S24 is the flux density measured within the 24-flm MIPS 
band, and we explicitly denote SFRs calculated from this band as 
1*24 ■ We have not included any systematic contribution to the errors 
from uncertainties in this conversion between flux density and SFR. 

As with other methods of estimating SFR, there is an uncer- 
tainty due to the conversion between an observable parameter that 
comes principally from high-mass stars (> 5 Mr.,, which carry out 
the majority of the dust heating), and a SFR which includes the 
rate of formation of much lower mass stars (down to 0.1 M©). The 
Initial Mass Function (IMF) used to convert between the two is typ- 
ically ass umed to be a s ingle power-law slope of —1.35 from 0.1 to 
100 M Q dSalpetedll955l) . which can lead to an over-estimate of the 
number of low-mass stars which are forming. [Rieke et aT I d2009h 
use a modified IMF with a shallower slope at low mass (—0.3 fr om 
0.08 to 0.5 M® and -1.3 from 0.5 to 100 Mm ; lKroupall200ll) in 
order to calculate their SFR estimates, which results in a SFR that 
is a factor of 0.66 times the value which would be calculated from 
a Salpeter IMF containing the same mass of high-mass s tars. This 
factor has already been applied to the Rie ke et al] d2009h relation- 
ship in equation {2j, and we adopt it throughout the re mainder of 
this w ork where appropriate - in particular, to modify the Kennicutt 
d 19981) SFR relationships, which were originally calculated using 
a Salpeter IMF. Error estimates on our results do not include any 
contribution from uncertainties in the IMF. However, all tracers of 
star formation that we consider in this work (Ha, UV, IR) are most 
sensitive to high-mass stars, and a change in assumed IMF at the 
low-mass end will make a very small difference to the overall con- 
clusions, as it would alter the normalisation factors in the estimates 
of total SFR approximately equally at all wavelengths. 

We indicate the SFRs that correspond to a given 24-flm flux 
density on Fig. [3] for reference. The diagonal lines indicate the 
Kennicutt ( 1998) relationship between Ha flux and SFR, assuming 
different levels of e xtinction for all ga laxies; as has been seen in 
previous works (e.g. lDopita et al.ll2002r) . calculated extinction val- 
ues vary significantly from galaxy to galaxy. The stacked data has 
a typical extinction of between 1 and 2 mag for all flux bins, and 
a visible increase for higher SFRs. Note that the individual galaxy 
detections are biased against low extinction (and therefore low 24- 
flm SFR) sources since only Ha emitters with 24-flm counterparts 
are plotted - the stacked data has no such bias, as it uses both non- 
detections and sources individually detected at 24 flm. 



4.3 Variation in Ha extinction with observed Ha luminosity 

We can calculate the typical Ha extinction for a bin of data, if 
we are willing to assume that 24-flm emission for an Hq-selected 
sampl e is a reasonable SFR indicator, a s was found by Rie ke et all 
(2009). We invert the iKennicuttl d 19981) relationship between SFR 
and Ha luminosity (modified by a factor of 0.66 to make the IMFs 
consistent; see Section l4~2t to obtain a typical intrinsic Ha lumi- 
nosity, Ln a i nt , for each bin, i.e. 
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As galaxies are all at a known redshift, these luminosities are con- 
verted to the observed line fluxes by assuming that a single value 
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Figure 4. Variation in Ha extinction for the Ha emitters (black circles with error bars), after binning them by their observed Ha flux, (a) Variation in extinction 
with observed Ha luminosity (corrected for N II contamination but not for dust attenuation). The solid line shows the predicted relationship from equation (5), 
with the approximation from equation HO) shown by the dashed line. The best fit to the data, from equation |6), is indicated by the dotted line, (b) Variation 
in extinction with 24 - (im SFR. The solid line shows the predicted relationship from equation Q, using the fit be tween the Balmer dec rement and SFR from 
Hopkins et al. 12001), while the dashed and dash-dot lines are the equivalent relationships, using the fits from Sul livan et alj 4200lh andlPerez-Gonzalez et al. 
12003) respectively - see Section |431 for further details. The best fit to the data, using the slope of equation (7) and B = 5.8, is shown by the dotted line. The 
hypothesis of a constant extinction for all Ha luminosities or SFRs is excluded at the 3o~ and 4o~ level respectively. 



of extinction Ah« can be applied to galaxies in each bin, i.e. 



Ana = 2.5 log I0 



J Ha,int 
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where is the luminosity-distance at z '- 
equations, we obtain 



- 0.845. Combining these 
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We calculate an extinction from each set of binned data using equa- 
tion 0, and the results are shown on Fig.[4fa). There is an indica- 
tion of an increase in Ah« Wlin observed Ha luminosity; the best 
fit to the data is 

A H « = (- 19.46 ±6.77) + (0.50 ±0.16) log 10 ( ^^f] , (6) 

V erg s 1 / 

i.e. a 3o" detection of a variation in extinction with observed Ha 
luminosity. However, if there is a dependence on luminosity (and 
by implication, on SFR activity) then relating extinction to ob- 
served Ha luminosity may obscure some of the effects. In or- 
der to test this, we show the dependence of extinction on SFR in 
Fig. |4fb). The dependence of extinction on SFR is more promi- 
nent, with a best fit to the data of A H « = (0.73 ± 0. 15) + (0.44 ± 
0.11) log 10 (1*24 /^0 yr -1 )- i-e- a 4rj detection of SFR-dependent 

extinction. 

P revious studies (e.g.lW ang & Heckman 1996; Hopkins et al. 
200 ll; ISullivan et al.l)200ll; iDopita et alj l2002r. iBerta et alJ^Ool 
Perez-Gonzalez et al.ll2003l ; iBuat et alj|2005l ; ISchmitt et alj|20061 ; 
Caputi et alj 120081) have found a correlation between extinc- 
tion and either luminosity or SFR, with the implication be- 
ing that mo re actively star - formi ng regions within a galaxy are 
also dustier. iHopkins et al. (2001) derived a SFR-dependent red- 
dening relationship, based upon the observed Balmer decrement 
and far-IR luminosity of a sample of local star-forming galaxies 
dWang & Heckman 1996). We repeat the derivation from Section 3 
of lHopkins et alj ( 120011) in order to obtain a general expression for 
their relationship between SFR and extinction, 



'Ha 1 



2-5£ H a 



lo gl() 



0.797 log IO (Tj>P)+3.83 



(7) 



where is related to the E(B — V) colour excess by = 
Ax/E(B — V), and kna and k^p are the values of k^ at the wave- 
lengths of Ha and H/3 radiation respectively. *P is the total SFR 
of a galaxy in Mq yr -1 , and rj is a conversion factor appropriate 
to the chosen IMF (1 for a Salpeter IMF, and 0.66 for the Kroupa 
IMF used in this work). The factor of tj in equation cancels an 
equivalent factor of \jr\ during the calculation of from IR lu- 
minosity in Section |4~2l and is otherwise unimportant. The value 
of 2.86 represents the intrinsic Ha/H/3 line flux ratio (appropriate 
for Case B recombinat ion, a temperature of 10 4 K, and an electron 
density of 10 cm" 3 ; lBrockIehurstlll97lh . and the other numerical 
factor s come from the least-squares fit performed by Hopki mTet alj 
d200ll) to the observed Balmer decrement and the logarithmic far- 
IR luminosity of local star-forming galaxies. 

In order to test the form of the dust attenuation law, we assume 
that the numerical factors within the square brackets of equation 
are fixed and correct, and fit our data for the numerical prefactor, 
B = 2.5kna /(knp ~ ^Ha)> with errors coming from the typical un- 
certainty inAjr a . We obtain B = 5.8 ± 1.0, or alternately 



K np 

fctta 



= 1 



2.5 



= 1.43±0.07. 



(8) 



This i s consistent with the value obtained from the ICalzetti et alj 
(2000) dust attenuation law, which predicts a value of B = 6.54, and 
gives £ H| g /kna = 1-38. We show this result on Fig.|4[b) to demon- 
strate the good agreement with our data. Values of B obtained from 
the i ndividual fields are also c onsistent with this value (see Tabled- 
iBrinchmann et al. U2004h have tested the two assumptions that 
we have made in the derivation of equation 10, namely that the 
intrinsic Ha/H/3 ratio is always equal to 2.86, and that the con- 
version between intrinsic Ha luminosity and SFR that we use is 
appropriate for all galaxies. They conclude that while both assump- 
tions are slightly flawed - the Case B ratio depends slightly upon 
temperature, while the most metal-rich galaxies produce a slightly 
lower Ha luminosity for the same SFR - these variations are small 
and approximately cancel each other out; the combination of the 
two assumptions is therefore a fairly good approximation to make. 
Variations in the intrinsic Case B ratio for different t emperatures 
(between a ratio of 2.72 and 3.03; iBrocklehursil 197 lh could only 
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affect our expected value of B = 6.54 by up to ±0.16, which is 
much smaller than t he stat istica l errors in this study. 

ISullivan et all ( 1200 lh and IPerez-Gonzalez et alj d2003h have 
also produced least-squares fits to the correlation between the 
Balmer decrement and Ha-derived SFR from independent sets 
of galaxies, leading to slightly different numerical factors within 
the square brackets of equation . We show the equivalent ver- 
sions of equation 0, calculated from these fits, on Fig. |4fb) for 
comparison with th e IHopkins et al.l d200lh result. As noted by 
ISullivan et alj (l200ll). the slope of their correlation is very similar 
to that of IHopkins et "all feOOlh . although the indep endent estima- 
tion o f SFR from IR observations (carried out by Hopki nTet alj 
l200lL and in this work), is preferable to using Ha-derived val- 
ues, providing that AG N contamination has been removed. We find 
that the slope of the lSullivan et al relationship agrees with 

our data, but that the overall normal isation is too high. In con- 
trast, the IPerez-Gonzalez et al.1 {2003) correlation does not agree 
with our results; they indicate that if highly-obscured objects are 
re moved from their sam ple, they find a correlation similar to that 
of lSullivan et"all d200ll). 

We note that 1 Afonso et al. I J2003b have found a significantly 
steeper relationship between SFR and Ha extinctio n from a radio- 
selecte d sample than tha t obtained in th i s work or bv lHopkins et al.1 
d200ll) . However, while IHopkins et al, I booih found a reasonable 
correlation b etween the Balmer d ecrement and IR luminosity of 
their sample. lAfonso et"al] (2003) did not find a good correlation 
between the Balmer decrement and radio-derived SFR of their 
sources, which may make their derive d relationship less reliable. 

Following Hopk ins et al. I booih we combine equations {3} 
and 10 to obtain a relationship between the extinction of a galaxy 
and its observed Ha luminosity, 



A Ha = 12.8 xlO 4 ""- 2.5 log I0 



( Lh<x, 



7.9 x 10~ 42 I - 12.0, 

V erg s 1 J 

(9) 

where the observed Ha luminosity Lna.obs is me value calculated 
from the observed Ha line flux, after corrections for N II flux con- 
tamination, but before any dust attenuation correction has been ap- 
plied to the data. Equation ((9} cannot be solved analytically; how- 
ever, over the luminosity range covered by sources in this study 
(1.6 xlO 41 <L H «, obs 
by the linear fit 



< Lh« „k < 4 x 10 42 erg s '), we can approximate it 



Aua = -21.963 + 0.562 log I0 



erg s 



(10) 



This is consistent with our fit to the data from equation J6), as 
shown on Fig. EJa), and is within 5 per cent of the numerical 
solution over the observed luminosity range tested in this work. 
The approximate and numerical solutions are consistent to within 



10 per cent over the luminosity range of 10 — 10 erg s 

4.4 Dependence of Ha extinction on galaxy properties 

We have shown that equation is appropriate for our full sam- 
ple of Ha emitters, and accurately describes the variation in dust 
extinction for sources with different SFRs. In order to test whether 
there is any variation in dust attenuation between different types of 
galaxies, we split the sample into subsets based upon other prop- 
erties of the galaxies. We stack each subset into different Ha flux 
bins in order to remove the dependence of extinction on SFR, as the 
mean SFR of different galaxy types can vary strongly. We then re- 
peat the analysis of Section l4~3l in order to verify that equation 



Table 2. The calculated values of normalisation B and error c% for equa- 
tion {7}, broken down by the source selection criteria. The number of galax- 
ies that satisfied a particular selection criteria, N, is listed for reference - if 
N was not large enough for meaningful statistics, no normalisation is given. 



Sample 


N 


B 




All Ha emitters 


618 


5.8 


1.0 


COSMOS 


422 


6.2 


0.8 


UDS 


196 


4.9 


1.5 


COSMOS: 








disk / spiral morphology 


347 


6.9 


0.7 


irregular morphology 


57 


6.3 


2.4 


early-type morphology 


12 






disk / spiral, disk-dominated 


117 


6.5 


1.1 


disk / spiral, intermediate 


155 


7.3 


0.7 


disk / spiral, bulge-dominated 


21 






merging 


92 


6.0 


1.9 


not merging 


291 


6.9 


0.7 


low-density 


207 


6.8 


0.8 


high-density 


207 


6.3 


0.9 


9.5 < log 10 (stellar mass/M ) < 9.8 


67 


1.6 


1.6 


9.8 < logio(stellar mass/M ) < 10.1 


141 


4.6 


0.9 


10.1 < logi (stellar mass/M Q ) < 10.4 


180 


6.6 


0.6 


10.4 < logi (stellar mass/M ? ) < 10.7 


135 


7.8 


0.7 


10.7 < log| O (stellarmass/M ) < 11.0 


57 


6.8 


2.7 


Predicted value from Calzetti et al. ("2000) 




6.54 





accurately describes the variations in the data, and to calculate a 
value of B for the subset. A summary of the relationships derived 
for all subsets of galaxies can be found in Table|2] 



4.4.1 Morphological type 

A large Hubble Space Telescope survey of the full 1.8 deg 2 COS- 
MOS field has been completed dScoville et ai] |2007b), producing 
images with ~ 0. 1 arcsec resolution. Morphological classifications 
have been performed on ~ 40 per cent of the field dScarlata et al.1 
2007), and the majority of the Ha emitters have b een assigned 
a sourc e type (disk/spiral, irregular, early-type) by ISobral et al.1 
(2009a), based upon this morphological classification, visual in- 
spection, and colours taken from the Subaru optical data. There are 
too few early-type galaxies in the sample for a meaningful compari- 
son, but the disk/sp i ral and irregular galaxies are shown in Fig.[5ja). 

IScarlata et al] d2007l) have further classified the galaxies with 
a disk/spiral morphological type by their bulge-to-disk ratio. The 
number of bulge-dominated galaxies is too small for a meaning- 
ful comparison, but no significant difference is seen between the 
intermediate and disk-dominated galaxies in Fig.|5lb). 



4.4.2 Merger activity 

ISobral et al] d2009al) have visually inspected the galaxies and clas- 
sified them according to whether they are undergoing mergers (93 
galaxies), or definitely not merging (292 galaxies). The remaining 
galaxies were potential mergers; we have not included these in the 
analysis as it is unclear which class they should belong to. The 
merging galaxies and non-merging galaxies are shown on Fig.[5]c). 
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Figure 5. The relationship between Ha extinction and SFR, for various subsets of the COSMOS data - the meaning of each classification can be found in 
Section |4~4l The diagonal lines show the predicted relationship between An a and SFR from equation using B = 6.54, and are not fitted to the data. Values 
of B obtained through fitting equation (7) can De found in Tableff] along with errors for each fit. 



4.4.3 Environment 

The angular distances to all nearby sources at comparable redshift 
have been calculated for each Ha emitter, and converted into a lin- 
ear distance to the lOth-nearest neighbour, r\$ (D. Sobral et al., in 
preparation). From these distances, a local projected surface den- 
sity can be calculated using £ = 10/7Tr 2 () . We split the sample in 
half according to the local density; as shown in Fig.[5jd), no signif- 
icant difference is seen between the two subsets. 

4.4.4 Stellar mass 

Stellar mass estimates are av ailable for all COSMOS and UDS 
sources dMobasher et~ai1l2007l . M. Cirasuolo et al., in preparation). 
We split our sample in half by stellar mass, as shown in Fig.[6ja). 
Two results are apparent from this plot - bins with lower stellar 
mass also have a lower median SFR, and the lowest stellar mass 
bin appears to have a decreased normalisation compared with the 
higher mass bin. 

Previous studies have sho wn that the stellar mass and SFR 
of a galaxy are correlated (e.g. iBrinchmann et alj|2004T) . and the 
first of our results is consistent with these studies. In order to in- 
vestigate the second result further, we divide our sample into five 
stellar mass bins (and three Ha flux bins), calculate normalisations 
for each sub-sample, and use equation 10 to predict the extinc- 
tion that would be obtained for a galaxy of this stellar mass, and 
a SFR of 10 Mr. } yr -1 . These extinctions are shown in Fig.[6lb). 
We find that galaxies with stellar mass above ~ 10 10 ' 2 Mq have 
an approximately constant extinction / SFR relationship (i.e. for a 
given SFR, the predicted extinction would be the same), but that 
below this stellar mass, the typical extinction for a galaxy with 
given SFR decreases significantly. This may be a result of incom- 
pleteness in our sample, based upon the stellar mass distribution of 
our Ha emitters (which we find to peak at ~ 10 10,2 Mq). We do 



however note that Brinchm ann et alj ( 120041) have detected a depen- 
dence of Ha extinction on stellar mass for galaxies within SDSS, 
which they attribute to a correlation between mass and metallic- 
ity for star-forming galaxies. Their results do not take into account 
the dependence of extinction on SFR which we have confirmed in 
this work, and so may simply be a consequence of a correlation 
between SFR and stellar mass for star-forming galaxies. We defer 
further investigation of this relationship to a separate work (Garn et 
al. in preparation). 



4.4.5 Summary 

With the exception of the potential stellar mass dependence, we find 
no statistically significant difference between the normalisations of 
the relationship between extinction and SFR for any of these sub- 
sets of galaxies, as shown in Tableland Fig. [5] This may simply 
be a consequence of the limited number of galaxies available - af- 
ter splitting the COSMOS galaxies into two or three classifications, 
and a further three Ha flux bins, there are only typically a few tens 
of galaxies contained within each bin, and any difference would 
have to be very large in order to be detectable. Note that plots of 
extinction against observed luminosity (not shown) also indicate 
that there is no significant difference between the various subsets, 
but we have chosen to illustrate this fact by relating the extinction 
to 24-jUm-estimated SFR rather than to o bserved Ha luminosit y, so 
that a direct comparison to the results of Hopki ns"et al] d200lh can 
be made. 

We conclude that within the uncertainties quoted in Table [2] 
the relationship between extinction and SFR that was presented in 
equation (|7) is appropriate for use in galaxies with a wide range of 
morphologies and properties. 
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Figure 6. The relationship between Ha extinction and SFR for data binned by stellar mass, (a) Data separated into two stellar mass bins, (b) The extinction 
that would be predicted for a galaxy of a given stellar mass, and a SFR of VP = 10 Mq yr - 1 , with the data separated into five mass bins. 



Table 3. The rest-frame wavelength, zero-point correction and median 
galactic extinction for each of the three bands used to calculate the UV 
extinction, along with the calculated conversion factor y between nebular 
and stellar reddening for each band from equation )13t . using the stacked 
data (Fig. [7), and from individual sources (Fig. [8). 





u* 


Bj 


8 + 


Rest-frame X (A) 


2058.5 


2417.2 


2590.6 


Zero-point correction 


-0.084 


0.189 


-0.090 


Galactic extinction 


-0.086 


-0.074 


-0.069 


7 (stacked data) 

7 (individual sources) 


0.50 ±0.07 
0.49 ±0.14 


0.50 ±0.07 
0.49 ±0.14 


0.48 ±0.07 
0.47 ±0.14 



4.5 UV dust extinction 

The UV luminosity of a galaxy is also produced directly from 
young, high-mass stars, and will be attenuated by interstellar dust. 
We estimate the dust extinction in the rest-frame UV, A\jy, in a 
similar manner to Section l4~3l As before, we assume that 1*24 is a 
good proxy for the true tota l SFR, and use an IMF-modified ver- 
sion of the[Kennicutt ill 9981) calibration to calculate a typical UV 
luminosity, L v , for each bin, through 



Ly 



WHz- 



:4.71 x 10 



*24 



(ID 



V M Q yr -1 

Ajjv is then estimated through comparison of L v and an observed 
flux density, S v , using 



A uv = 2.5 log! 



Ly(l- 



(12) 



where the 1 + z factor accounts for bandwidth stretching. 

We have three photometric bands that trace far-UV flux den- 
sity (u*, Bj and g + ), with the rest-frame wavelengths that they 
probe listed in Table [3] We only consider COSMOS galaxies in 
this section, as there are currently no accurate photometric correc- 
tions available for the UDS data. We applie d the photometric ze ro- 
point corrections that are recommended by Cap ak et all d2007t) to 
the photometry of each source, along with individual corrections 
for galactic extinction (with the median correction factors listed in 
Table[3]l. An appropriate ape rture correction (which applies to each 
band; see Capa k et al.||2007t) was taken from the PSF-matched cat- 
alogue for each source, with a median correction of —0.275 mag. 
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Figure 7. Variation in the median UV extinction with 24-/im SFR, for the 
HiZELS COSMOS emitters after selection by their Ha flux. The UV ex- 
tinction has been calculated from photometry in the u* band, B j band (off- 
set by +1 mag) and g + band (offset b y +2 mag). The solid lines are fits to 
equation U3K using the Calzetti et al. (2000) dust attenuation law to calcu- 
late k\jy, and with values of y as listed in TablefJ] 



Fig. [7] shows the measured dependence of UV extinction on SFR 
for the three bands - as before, there is a correlation between the 
SFR and UV extin ction. 

ICalzettildl997h found that the dust obscuration of stellar con- 
tinuum emission between the wavelengths of the Ho: and H/3 lines 
is less than would be expected from an extrapolation of the obscu- 
ration of the line emission. They found that the optical depth of 
the dust that was obscuring the continuum was about 60 per cent 
of the depth of the dust that obscures the line emission. The infer- 
ence drawn from this was that massive stars (producing the Ha/H/3 
emission) form in dusty environments, while the bulk of the red 
stellar continuum was a result of much less dust-obscured stars 
(presumably older, and lower mass), at a different location in the 
galaxy. They suggest using E(B - V) sW = (0.44 ± 0.03) x E(B - 
V) gas to calculate the reddening of the continuum emission. 

The difference of environment should not be true for Ha and 
UV continuum radiation, which are both predominantly produced 
by the same population of young stars; it might therefore be ex- 
pected that the correction factor (which we denote as y) should 
converge to 1 for shorter wavelength continuum radiation, where 
the obscuring dust will become more like the dust that obscures the 
Ha emission. Not e, however, that the v alue of y is closely linked 
to the shape of the ICalzetti et al.l ( feOOOl) dust law, which has been 
calculated through comparisons of the SEDs of attenuated galaxies 



12 



T. Gam et al. 



0.2 



0.1 




0.5 



7 



Figure 8. The distribution of y for COSMOS Ha emitters with individual 
detections at 24-/im, calculated from the u*,Bj and g + bands. The Gaussian 
is a fit to the M-band data, with y = 0.50 and a = 0. 14. 



and their unattenuated counterparts - different wavelength regions 
of the dust attenuation law relate to different regions of a galaxy 
(i.e. the UV end is physically related to the dusty star-forming re- 
gions, while the long wavelength end will be weighted towards the 
old stellar population residing in less dusty regions). As such, there 
is a degeneracy between the wavelength dependence of y and kx 
that does not relate to physical properties of dust within galaxies. 
We therefore believe that while the value of y is useful in calcu- 
lating UV attenuation from the Balmer decrement and vice versa, 
it does not have a direct physical relevance, as it is closely linked 
with the assumed form of the dust law. If UV extincti on at z = 0.84 
is the same as that seen locally, and described by the Calz etti et al.l 
( 2000) dust law, we would expect to find y ~ 0.44 for UV emission 
as well as for longer wavelength continuum emission. 

We allow the conversion between nebular and stellar redden- 
ing, y, to be a free parameter, and adjust equation 10 by this factor 
to obtain 



2.5y£ uv 



log] 



0.797 log 10 (?7'P) + 3.83 



2.86 



(13) 



where k\jy is me dust attenuation at the rest-frame wavelength of 
interest. We fit for yby assuming that the lCalzetti et al.l OoOOh dust 
attenuation law accurately estimates k\jy (which we verify through 
a comparison of the results from different bands), and show these 
fits on Fig. [7] The values of y are gi ven in Tablejl]- all are larger 
than 0.44, but are within ~ la of the lCalzettil ( 1997) result. 

We calculate the Ha and UV extinction for each of the 134 
COSMOS Ha emitters with a 24-/xm counterpart in the MIPS cat- 
alogue, that were not flagged as AGN. A strong correlation is seen 
between Au a and eacn of the Auv measurements (with p = 0.70), 
and from these extinctions we calculate y on a source-by-source 
basis. The distribution of yis shown in Fig. [8] and the results given 
in Table[3]- in each case the data is consistent with a Gaussian cen- 
tred on 0.5, with a = 0.14. A two-sample Kolmogorov-Smirnov 
test confirms that no significant difference is seen between the 
distributions of y in the three bands. We therefore conclude that 
y = 0.50 ±0.14, i.e. the dust attenuation of UV emission is ap- 
proximately half of that which would be pre dicted from a simple 
extrapolation of nebular reddening using the lCalzettiet al] ( bOOfj) 
dust attenuation law. Th is is in agreement w ith the low-redshift re- 
sult, and implies that the lCalzetti et al. (2000) law is broadly appli- 
cable to our galaxies. 

The principal qualitative difference between the lCalzetti et al.l 

(2000) dust attenuation law, derived from local starburst galaxies, 
and the individual extinction laws measured from observations of 



stars within the Milky Way or the Large Magellanic Cloud is the 
absence of a 'bump' at around 2175 A in the rest-frame UV. There 
is a conceptual difference between dust 'extinction' by an interven- 
ing dust cloud, and dust 'attenuation' from an exten ded source by a 
distribution of dust that is m ixed with the stars (see Calzetti 1 19971 
for more details). However, iGordon. Calzetti & Wittl h997l) have 
shown that is it not possible to explain the absence of a 2175 A 
dust feature purely through geometrical effects, implying that the 
lack of this feature is intrinsic to the properties of dust grains within 
local starburst galaxies. 

There are indications that a few high-redshift sources may 
contain this dust f e ature (e.g. ISolorzano-Inarrea et alj [20041: 
I Srianand et al] 120081 : lElfasdottir et al] 120091 : iNoterdaeme et al] 
2009). A 2175 A bump would show up prominently in our data, 
via a greater extinction in the «*-band data than in the Bj and g + 
bands. As no such difference has been found, we conclude that the 
dust attenuation law of high-redshift star-forming galaxies does not 
require the addition of a UV bump. 



5 DISCUSSION 

5.1 Variation in Ha dust extinction with SFR 

Ha emission is an excellent tracer to use for cosmic SFR studies 
as it is sensitive, can be easily studied over large volumes, and has 
a clean selection function. However, it is necessary to correct the 
observed Ha luminosity for the effects of dust extinction before 
the true SFR can be obtained. 

We have shown in Section|4]that the empirical relationship be- 
tween the SFR of local galaxies and their Ha extinction is also ap- 
propriate for galaxies at z = 0.84. We have further split our sample 
of galaxies by their morphological type, bulge-to-disk ratio, merger 
activity and the local density of their environments. No statistically 
significant difference was seen in the normalisations of the SFR / 
Ha extinction relationship, implying that this is a general property 
of galaxies which applies at least over the range < z < 0.84, to 
spiral and irregular galaxies, disk-dominated galaxies and galaxies 
with a moderate bulge component, to galaxies that are and are not 
undergoing merger activity, and galaxies in high and low density 
environments. 

The increase in extinction with higher SFR can be explained 
through the fact that star formation takes place within dusty molec- 
ular clouds, and the SF R is related to the gas surface density (the 
Schmidt-Kennicutt law: ISchmidtll959l : lKennicuttll989l) - more ac- 
tively star-forming regions are likely to have dustier environments, 
where the Ha luminosity produced by young stars will suffer a 
greater amount of attenuation. 

Equation Q allows us to estimate the typical extinction that 
is appropriate for galaxies with different SFRs; the value of Au a = 
1 mag - commonly taken to be the amount of dust extinction to cor- 
rect for - is appropriate for galaxies with a SFR of ~ 3 Ma> yr , 
comparable to the typical rates of star formation in the local Uni- 
verse. The typical extinction that is appropriate for galaxies with a 
SFR of 10 and 100 M yr~' is ~ 1.3 and ~ 1.7 mag, although ex- 
tinctions calculate d for individual gala xies may vary considerably 
from these values. iDopita et al] <2002h find a range of Ha extinc- 
tions between ~ 1 and 3 for galaxies with SFR of ~ 100 M@ yr , 
and Fig.|3]shows a similar spread of IR-derived Ha extinctions for 
individual sources. 
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5.2 Variation in Ha dust extinction between < z < 0.84 

Several works (e.g. Ilresse etaDl200l fvillar et al.ll2008h find that 
the amount of dust e xtinction is greater at high redshift than in the 
local Universe, with IVillar etal.U2008h obtaining a mean extinc- 
tion of An a = 1 .48 mag from a sample of Ha emitters at z = 0.84. 
We also find this increase, with a typical extinction value within 
1<T of 1.5 mag; however, our results from this section and Sec- 
tion [43] demonstrate that this merely reflects the greater number 
of high-SFR sources at higher redshift, rather than an intrinsic in- 
crease in extinction for similar galaxies at the two redshifts. We 
have excluded the possibility of there being no dependence of Ha 
extinction on SFR at the 4<J level, a concl usion that is in agre ement 
with the results of IChoi et all fe006h and lCaputi etal] fe008h . who 
find that variations in V-band extinction for different galaxies are 
correlated with luminosity, rather than with redshift. 

The typical extinction that is measured from a survey that de- 
tects Ha emitters at a specific redshift will be dependent on the 
detection threshold - counter-intuitively, a deeper survey will lead 
to a lower measurement of the average extinction, since from Fig.[4] 
a lower luminosity threshold also implies a lower extinction thresh- 
old. We use this result to caution against drawing conclusions about 
the typical dust extinction of galaxies at a given redshift, without 
taking into account any SFR dependency. 



5.3 The effect of variable extinction on the Ha luminosity 
function and star formation rate density at z = 0.84 

ISobral et alJ l[2009a) presented an extinction-corrected Ha lumi- 
nosity function at z = 0.84 from the HiZELS data used in this work. 
A single extinction of An a = 1 mag was used to correct all sources, 
and a Schechter function was fitted to the data, yielding values of 
L*, <j>* and a which are presented in column 1 of Table|4] 

We recalculate the Ha luminosity function at z = 0.84 from all 
Ha emitters, applying the extinction correction from equation J 1 0b . 
and obtain new values of L* , tf>* and a, listed in column 2 of Ta- 
ble [4] These results remain consistent with s trong evolution in L * 
and <j>* over the range < z < 1, as found in ISobraletal.l (l2009a). 
The most significant change is an increase in the characteristic 
turnover luminosity, L*\ the increase in the typical Ha extinction 
leads to a greater number density of bright emitters, compared with 
taki ng A^ a = 1 mag for all sources. 

ISobraletall l 2009a) calculated the average star formation rate 
density (SFRD) at z = 0.84, by integrating the luminosity function 
and applying an average 15 per cent correction for AGN contam- 
ination. If we use a variable Ha extinction correction, and simply 
apply a 15 per cent AGN contamination correction, we obtain a 
significantly larger value for the SFRD than that found from us- 
ing 1 mag of extinction (column 2 of Table 0. However, we can 
improve upon this estimate by recalculating the Ha luminosity 
function, excluding those sources which were classified as potential 
AGN (column 3 of Table EJ. There is a flattening seen towards the 
faint end of the luminosity function, with a consistent with a value 
of —1.4. We integrate this luminosity function down to the limit 
of the HiZELS survey, L Ha = 10 4L5 erg s" 1 , and obtain a SFRD 
of 0.15±0. 01 Mp) yr -1 Mpc~ 3 , essentially identical to the value 
obtained bv rSobral et all b009al) of 0. 15 ± 0.02 Mr, yr -1 Mpc~ 3 . 
By integrating the entire luminosity function, we obtain a SFRD 
of 0.32 ±0. 1 2 M B yr -1 Mpc~ 3 , compared with the value from 
ISobral et al.l J2009j) of 0.37 ±0.18 M G yr" 1 Mpc~ 3 . The com- 
bination of an improved AGN identification method and variable 



Table 4. Schechter function parameters L*, 0* and a for the z = 0.84 Ha 
luminosity function, along with the average SFRD (in Mj, yr~' Mpc~ 3 ) 
obtained by integrating the luminosity function down to the HiZELS sur- 
vey limit of Luc = 10 41 ' 5 erg s~', or across the total luminosity range. The 
Ha exti nction assumed whe n calculating the luminosity function is either 
'1 mag ' JSobral et al]2009 j) or 'Variable' (the results from this work), and 
the AGN correction applied when calculating the SFRD either uses all Ha 
emitters to calculate the luminosity function and then assumes a 15 per cent 
contamination fraction, or only uses those emitters which were not identi- 
fied as AGN. 



Corrections: 








Extinction 


1 mag 


Variable 


Variable 


AGN 


15 per cent 


15 per cent 


AGN excluded 


Parameter: 








logio(L*/ er g s ~') 


42.26 ±0.05 


42.50 ±0.04 


42.44 ±0.05 


logio(07Mpc- 3 ) 


-1.92±0.10 


-2.05 ±0.09 


-1.99±0.10 


a 


-1.65 ±0.15 


-1.60±0.10 


-1.40 ±0.20 


SFR density: 








> 10 4L5 erg s" 1 


0.15 ±0.02 


0.24 ±0.05 


0.15±0.01 


Total 


0.37±0.18 


0.41 ±0.16 


0.32±0.12 



dust attenuation only has a moderate effect on the calculated SFRD 
at z = 0.84, although may be more important at higher redshift. 



6 CONCLUSIONS 

We have used a well-defined sample of Ha-emitting galaxies at z = 
0.84 to study the link between star formation and dust extinction. 
The HiZELS sample consisted of 694 Ha-emitting galaxies within 
the redshift range z = 0.845 ±0.021 in the well-studied COSMOS 
and UDS fields, with additional information available in multiple 
UV, optical and IR bands, along with ancillary radio and X-ray data. 

In order to focus on star-forming galaxies, it was necessary 
to remove AGN contaminants from our sample. We performed this 
removal through a number of multi-wavelength classification meth- 
ods; emission-line diagnostics, deviations away from the IR / radio 
correlation, X-ray counterparts, mid-IR and mid- to far-IR shapes 
and full template fitting were all used to identify potential AGN. We 
estimate an AGN contamination rate of between 5 and 1 1 per cent 
for Ha-selected emission-line sources at z = 0.84. No single di- 
agnostic identified more than half of the AGN contaminants, with 
SED template fitting being the most successful; we suggest that us- 
ing a combination of template-fitting and colour-colour diagnostics 
will maximise the success rate of AGN identification in Ha sur- 
veys. 

In order to constrain the average IR luminosity of Ha emitters 
at z = 0.84, we have binned our sample by observed Ha flux, and 
used the technique of source stacking to calculate typical 24-flm 
flux densities for the Ha emitters. Several previous studies have 
found that dust extinction correlates with either luminosity or SFR, 
and we looked for a relationship between the SFR of a galaxy and 
its dust attenuation. We found that individual extinctions vary sig- 
nificantly between galaxies, but confirm that this correlation holds 
for the average properties of our sample; galaxies with higher SFR 
also undergo a greater amount of extinction. 

We have presented an approximate method of calculating the 
typical Ha extinction of a galaxy, based solely upon the observed 
Ha luminosity; A Ha = -21.963 ± 0.562 log 1() (L Ha>obs /erg s~'), 



14 T. Garn et al. 



which is within 10 per cent of the full numerical solution over the 
luminosity range 10 41 — 10 44 erg s _1 . This is found to describe 
the properties of galaxies with disk, spiral and irregular morpholo- 
gies, a range of bulge-to-disk ratios, galaxies in high- and low- 
density environments, and both merging and non-merging galax- 
ies, although a potential deviation in this relationship may be seen 
for galaxies with low stellar masses. We emphasize that this is only 
an average relationship, and that the extinction measured within 
individual galaxies can vary away from this relationship by up to 
~ 2 mag. 

We have calculated the far-UV extinction of galaxies in a sim- 
ilar manner to the Ha extinction, and found a qualitatively similar 
correlation between this extinction and the SFR. We find that the 
UV dust obsc uration is about half of that expected from an extrap- 
olation of the lCalzetti et al.l j2OO0t) dust attenuation law. By calcu- 
lating UV extinctions from flux densities measured in three differ- 
ent bands, we show that, similar to the local Universe, there is no 
evidence for a 2175 A UV bump in the dust attenuation of high- 
redshift star-forming galaxies. 

The relationship between dust extinction and SFR at z = 0.84 
is the same as is seen at z ~ 0; we conclude that there is no evidence 
for a change in the dust properties of galaxies over this redshift 
range, and that any average increase in the extinction of galaxies at 
higher redshift is due to the overall increase in star formation activ- 
ity out to z ~ 2. Using our relationship between observed Ha lumi- 
nosity and extinction, we recalculate the Ha luminosity function 
for galaxies at z = 0.84. We find an increase in the characteristic 
luminosity compared with the value obtained by assuming a single 
value of extinction is applicable to all sources, implying a greater 
number density of bright Ha emitters than previous estimates at 
z = 0.84. 
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